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Abstract 

The reported CP asymmetry in D— > K + K~ /ir + Tr~ is argued to be too large to nat- 
urally fit the SM. If so, a new source of CP violation is implied in the AC = 1 sector 
with a milliweak strength. CP-odd interactions in the flavor-diagonal sector are 
strongly constrained by the EDMs placing severe limitations on the underlying the- 
ory. While the largest effects usually come from the New Physics energy scale, they 
are strongly model-dependent. Yet the interference of the CP-odd forces manifested 
in D decays with the conventional CP-even AC = 1 weak interaction generates at 
the charm scale a background level. It has been argued that the d n in the SM is 
largely generated via such an interference, with mild KM-specific additional sup- 
pression. The reported CP asymmetry is expected to generate d n of 30 to 100 times 
larger than in the SM, or even higher in certain model yet not quite natural ex- 
amples. In the SM the charm-induced loop-less |<i n | is expected around 10 -31 e-cm. 
On the technical side, we present a compact Ward-identity-based derivation of the 
induced scalar pion-nucleon coupling in the presence of the CP-odd interactions, 
which appears once the latter include the right-handed light quarks. 



1 Introduction 



Recent hints at possible direct CP violation in singly Cabibbo-suppressed D meson decays 
have caused some excitement, since they may be the first direct indication for physics 
beyond the Standard Model (SM). The LHCb collaboration observes a difference of time- 
integrated CP asymmetries [T] 

Aa CP = a CP (D° ->■ K + K~) - a CP (D° ->■ = -(0.82 ± 0.21 ± 0.11)%, (1) 

a result preliminarily confirmed by CDF [2J: 

Aa CP = -(0.62 ±0.21 ±0.10)%. (2) 

The CP asymmetry is defined according to 

(no^n r(£° -> f) - r(p° -» J) 

r(D°^/)±r(D ^/) 

Since the both final states are positive CP eigenstates in strong interactions, Aacp evi- 
dently roots in a direct CP asymmetry. 

As argued below, an effect of this magnitude solely in the framework of the SM may 
not be rigorously excluded, yet it would require a strong enhancement of certain decay 
matrix elements. Such a loophole definitely deserves further scrutiny. In this paper we 
analyze immediate consequences of the assumption that the reported effect is due to a 
new source of CP violation, beyond the CKM mechanism. 

In the present paper we focus on the impact of new CP-odd forces on the electric 
dipole moments, in particular on one of the neutron, d n . The observation at LHCb and 
CDF assumes CP violation in |AC| = 1 amplitudes. A more general model should embed 
this into a full flavor framework, hopefully highlighting a certain underlying symmetry. It 
would allow one to obtain predictions for other sensitive processes as well, in particular 
for B and K decays, where some tensions have also been noted in certain cases. 

We do not attempt to elaborate such a framework here and rather moderate the 
ambitions to noting that a flavor- diagonal CP violation of the size reported in the decays 
D° — y K + K~ and D° — > 7r + 7r~ is incompatible with the current limits on the electric 
dipole moment of the neutron. At the same time, we find that if the new CP-odd forces 
show up at low energies only in |AC| = 1 interactions, the neutron EDM is still well below 
the current limit, although it should be significantly enhanced, by more than an order of 
magnitude, compared to the SM. 

In the following section we briefly review the CP violation in D° -^-K + K~ and D° — y 
7r + 7r~ within the SM and introduce new |AC| = 1 interactions as a source of the enhanced 
CP violation. In Sect. 3 we examine d n in the SM and describe the loop- free mechanism 
to generate it at the charm scale. The elaborated estimates of the associated nucleon 
matrix elements indicate that it yields d n around 10~ 31 e-cm and may well constitute the 
principal contribution in the SM. The same analysis is then adapted to new BSM-mediated 
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CP-odd amplitudes to estimate the corresponding effect on the neutron EDM. Section 4 
summarizes the study. Appendix derives the induced CP-odd scalar pion-nucleon coupling 
by generalizing the Goldberger-Treiman relation and applying in QCD the current algebra 
technique. 



2 CP Violation in D° -> K + K~ and D° -> tt+tt 



2.1 Charm CP Violation in the Standard Model 

For the charm decays considered hereafter we have the effective weak interaction 



r - - Of. 

V2 



G, 



v cs v:-v cd v: d 



({cT^ulsT^-lcT^dT^d}) 



\ V cb V: b ([cT^sT^s] + [cT»u][dT,d] - 2[cT"u][&i>]) 



+ H.c. 



(4) 



= '-f: sin 9 C cos 9 C [oi - |r SM e n o 2 ] +H.c, 

V2 

[cT^iqT^q] = (c 7 ^(l- 75 )g) (^(l- 75 )«), T a 



7a(l-7s), 



where we have used the CKM unitarity, and color indices are assumed to be contracted 
within parentheses. The phase 7 is practically equal to the corresponding angle of the 
Unitarity Triangle, while 



?"SM 



viv ub 



v*v 



~ 7.5 x 10" 



(5) 



CP violation in the SM is quantified by the imaginary part of the invariant product of 
four CKM mixing elements describing the relative phase between the coefficients for 01 
and Oi-, 

A = Im V* s V us V cd V* d = Im {V^V cd )* {V* b V ud ) , (6) 

numerically A~3.3 • 10 -5 . 

The operators 01 and o 2 are a U-spin triplet and a U-spin singlet, respectively. Their 
interference induces the CP violation in the SM in the AC = 1 sector. In what follows 
we discard possible CP violation in D — D mixing since it drops out from the asymmetry 
difference Aacp- 

The decay amplitudes in D° — > f with / = / are given by (/ = K + K~ or f = n + n^) 



A(D° -> /) = -i-^ sin C cos 0c 
\/2 



m 



(/) 



(7) 



where m\ \ =(f\oi\D } are the reduced amplitudes, in general complex due to the strong 
interaction in the final state. The corresponding phases are generically referred to as 8\ ; 
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they are equal for the decays of D and D. Since the Oi amplitudes strongly dominate, 
fSM^lj the CP asymmetry takes a simple form 



a C p{D°-+f) = -r SM sin7 



(/) 



(/) 



sin 4?, r SM sin7~— — — -0.70-10 3 , 

sin Vc cos z Uc 



( f) (f) (f) 

where <5 21 =5\ — 5{ is the phase difference between the two hadronic matrix elements. 

The two final states K + K~ and 7r + 7r~ are components of the same U-spin triplet. 
Therefore in the SU(3) limit d^i 71 =^2\ K + 71 " would hold and 

|Aa CP | ~ 2\a CP (D ^K + K-)\ ~ 2|a CP ( J D°^7r + 7r-)|. (9) 

However, U-spin symmetry is significantly violated; one concludes from the decay rates 



T ^ pp ^ = Sr^^f-'S |mr ' 2 (10) 

that 

\mf +K ~ | ~ 0.456 GeV 3 , |mf^~ | ~ 0.252 GeV 3 . (11) 

We expect even larger potential SU(3) breaking in the phases of the amplitudes. 
The values in Eq. ffTTl) reasonably agree with the simplest factorization estimate 

m? +K ~ ^f^ K (M 2 K )f K (M 2 D -M 2 K ), mf*- ^-if^{Ml)U{Ml-Ml) (12) 

(the straightforward color renormalization factors have been omitted from the full expres- 
sion). It even yields the right scale for the S'f/(3)-breaking [3], although the literal ratio 
of the amplitudes tends to fall short of 1.81 in Eq. (II ip . 

The matrix elements of 02 determining the amplitudes m^ +i ^ and m^ 1 are not 
directly known. If estimated using factorization, one evidently obtains the values close to 
mi in Eqs. f fT2l . (II ip . with the additional minus sign for m^ 1 . However, the conventional 
factorization accounts only for the valence contributions. In a valence approximation, on 
the other hand, the same term in both o\ and 02 - with s-quarks for D° — > K + K~ and 
with (i-quarks for D° — >tt + tt~ , respectively - contribute. Consequently, no CP-asymmetry 
is generated in a valence approximation: the two strong amplitudes come from the same 
underlying operators and their strong phases coincide. 

Strictly speaking, any valence approximation should only be applied to the operators 
normalized at a low scale. The evolution of the operators 0\ and 02 above m& is identi- 
cal, yet generates additional terms for o 2 below it due to Penguin diagrams jl]. These 
in general have different strong phases. However, the new operators come with small 
loop-induced coefficients, while we do not expect their matrix elements to be enhanced. 
Therefore, we neglect these effects. 

As the starting point we assume that the magnitudes of mf +K and m^ 1 may be 
approximated using factorization, yet allow for arbitrary FSI phases relative to mi. This 
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amounts to having the ratio of the amplitudes in Eq. (jSJ) about unity. We then end up 
with 

\a CP (D°^K + K-)\ ^0.7 • 1(T 3 | sm5g +K ~\, \a CP {D ^0.7 ■ 10~ 3 | sin^" |, 

(13) 

and the sign of the two asymmetries may naturally be opposite. Therefore, the expected 
scale for |Aacp| in the SM is a few times 10 -4 up to 1.5 • 1CT 3 - provided the both FSI 
phase shifts are optimal. This is still about five times smaller than what is reported by 
LHCb. 

Accommodating the central value in Eq. (pQ) within the SM thus implies at least a 
five-fold enhancement of the U-spin singlet amplitude mediated by 02, or even a ten- fold 
if the asymmetry is dominated by one of the two modes and/or the strong phase shifts 
are not optimal. Moreover, this must happen for a non-valence part of the amplitude. 

Although at the moment the possibility of a sufficiently strong enhancement of the U- 
scalar amplitude in the D decay within conventional QCD dynamics cannot be rigorously 
ruled out, we view this possibility as contrived. The confirmation of the asymmetry Aacp 
at the currently observed level, in particular studying its share among the two channels, 
would be a strong evidence for the new CP-violating dynamics in the charm sector. At 
the same time, the strength of this conclusion crucially depends on the actual amount of 
the excess over our expectations. An eventual value around or somewhat below —0.3%, 
while still not smoothly accommodated in the SM, per se would make the case for new 
sources of CP-violation in D decays significantly weaker. 



2.2 Charm CP Violation through New Physics 

In what follows we adopt the assumption that the reported CP asymmetry roots in new 
CP-odd interactions. Within this hypothesis we will not attempt to stretch the uncer- 
tainties due to the QCD interaction to as strong extent and rather apply an educated 
judgment elaborated in weak decays so far; we then examine the consequences for the 
electric dipole moments. Neither we focus on the extreme values of parameters maximiz- 
ing the CP asymmetry. Consequently, we will gauge our expectations on an assumption 
that, speaking generally, the new source of CP violation produces an order of magnitude 
stronger CP-odd amplitude in D — » K + K~ or in D — > tt + tt~ decays than in the SM. 

Turning to NP, we make a relatively safe assumption that the New Physics-induced 
amplitude is small compared to the SM one this is obvious for the CP-odd NP part, 
and is applied also to its CP-even component. Then the asymmetry is given by the sum 
of the pure NP-induced asymmetry and the SM one, for either channels. Keeping in mind 
the conclusions of Sect. 12.11 we neglect the SM contribution altogether, and have 



a C p(D°^f) 



m 



(J) 



sin 8 



(/) 

NP J 



(14) 



where #np w|fp denotes the New Physics amplitude (in units of Gf sin 9c cos 9c /Vfy and 
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<5^p is its strong phase relative to of the SM0 

The value of the new couplings #np depends on the convention chosen to parameterize 
the BSM amplitudes (we tacitly anticipate using effective local operators to describe 
them). If we assume a 'natural' normalization of the operators where |?7ijyp| ~ Im^ | holds 
for the reduced amplitudes, we arrive at a ballpark estimate for the CP-odd coupling: 

|Im 5NP | ~ (2 5) • 10" 3 , (15) 

allowing for the generic unsuppressed strong phase differences as the educated guess about 
QCD dynamics in charm. The new CP-odd forces must therefore be of a 'milliweak' 
strength, according to the venerable terminology in CP violation. Their strength is in 
general about 10 times larger than what one estimates in the SM, cf. r^M sin 7 in Eq. ([8]). 
Specifically, to accommodate the 'direct' CP-asymmetry Eq. ([1]) one needs 

lmg NP m^~ lmg NP m^- + = -Aa CP 

10- 3 GeV 3 NP 10- 3 GeV 3 NP 8.2 • lCT 3 v ; 

Proceeding to the induced CP-odd effects in other observables requires specifying the 
nature of new interactions. Below the charm scale we have the realm of light hadrons 
including flavor-diagonal processes with stable hadrons and the decays of strange particles. 
CP-odd effects there are highly constrained; therefore we discard these, and relegate new 
sources to heavy particles. Their effect at low energies is described by local operators 
classified over the canonic dimension, with the lowest-dimension potentially dominating. 

A flavor-diagonal CP violation (for instance, the induced QCD #-term unless it is 
offset by a Peccei-Quinn-type mechanism) with a coupling of the size commensurate with 
Eq. ( fl5l) is by far excluded by electric dipole moments, in particular of the neutron. This 
likewise applies to the four-quark operators - they would generate d n in the ball park of 
10 -22 e-cm. Consequently, the flavor structure of the new CP-odd interaction must have 
vanishing flavor- diagonal components in the light sector. This property must replicate 
itself at the loop level, which strongly suggests it to apply to the heavy flavors alike. 
We do not analyze here the consequences of this requirement for various classes of the 
BSM models, but rather note that this may be a hint at an antisymmetric in generations 
structure of the underlying flavor dynamics. Yet even when postulating such a property, 
nontrivial constraints may follow at the loop level in view of the large gap between the 
scale of the potential effect on the EDMs and of their experimental limits. The loop- 
induced effects and the related renormalization of the effective low-energy operators may 
strongly depend on a particular class of models, see, e.g. Refs. [5j [6]; this lies outside the 
scope of our analysis. 

We therefore concentrate on the most direct consequences of the presence of new 
AC = 1 CP-odd amplitudes and describe them by the effective operators of dimension 5 
and 6. Most of them are four-quark operators. This appears representative enough. The 



1 We assume Iitl^np > and include the possible sign into Sj^ P . Let us also clarify that the phase 
convention required for CP conjugation is defined in such a way that the a\ amplitude in the SM is 
CP-even. 
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reason, as argued below, is that a significant - and probably the dominant - piece of d n in 
the Standard Model likewise originates from the same underlying effect: the interference 
of the CP-odd and CP-even weak |AC| = 1 amplitudes in the nucleon. The SM bears 
only mild additional model-specific suppressions; these may, or may not be vitiated by 
the new CP-odd |AC| = 1 interaction, depending on particular details. Consequently, we 
typically obtain a 30- to 100-fold enhancement of d n compared to the SM. 

The number of appropriate D = 6 operators (they can be both scalar or pseudoscalar) 
is quite large since they may differ in the chiral, color and light-flavor content. We first 
note that the scalar operators do not affect either of the D decays in question, yet they do 
generate d n . Therefore, d n could have been further enhanced if the scalar NP operators 
dominate. This possibility can be effectively eliminated experimentally by studying the 
similar CP-asymmetries including the parity-even final states in decays of D mesons, and 
we will not dwell on it any further. 

To substantiate the consideration we pick out ad hoc a few operators of interest: 

Oi = em c cia a pF al3 ^ 5 u , 2 = g s m c cia a pG aP ^ 5 u , 

3 = [cT^u]{\sT^s] + [d^d]), 4 = (c7 M (l+ 7B )u) (d^(l- l5 )d) (17) 



and put 



£ np = sin 9 C cos c ^2c k O k , (18) 



with Cfc dimensionless. The first two operators are the unique quark bilinears. Operator 
O3 has been picked since it evidently represents the CP-odd operator 02 of the SM - yet 
with an arbitrarily inflated coefficient. Consequently, we would roughly expect 

|Imc 3 | ~ 10- ir SM (19) 

if O3 is the only New-Physics source of CP violation. 

The operator O4 is an example with a different chiral content for both charm and light 
quarks and differs also in color and flavor. Operator O2, like O3 is a [/-spin singlet. For 
the sake of definiteness we assume in what follows that the direct CP asymmetry is largely 
seen in the 7r + 7r~ mode having a numerically smaller SM CP even amplitude, Eq. (II ip . 

The O4 matrix element can be estimated with simple factorization yielding 

1 2m 2 

(^-\O t \D) « -if,f°- (0)«S We (20) 

It would have shown a relative enhancement if charm mass scale were lower, while would 
have been suppressed for larger m c . For actual quark masses the corresponding factor is 
not too far from unity. This amplitude is color suppressed, therefore the factorization is 
not expected to be a good approximation - yet it makes explicit the expected qualitative 
features required for the scale estimates. 

The amplitudes for operators 0\ and O2 cannot be estimated by simple vacuum in- 
sertion. Keeping in mind that both are color-allowed we use instead a "rule of thumb" 
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for our estimates 



(n+n-\0'\D) T P T 



(tt+tt-\o\d) " V rp art ^ 

and set, as the reference, the operator O to be the 'valence' part of o%, (cT^d) (dT^u) 
(in fact, its P-odd part only). In other words, the fraction of the decay events into the 
exclusive tttc final state is assumed the same in the decays mediated by O and O' . Since 
charm mass lies in the intermediate domain, there must be no large kinematic factors 
floating around. 

For O2, the total decay width mediated by m c ug s ia^ v G^ v ^c is 



:N 2 -l 

jG 



a.m!^, (22) 



and the resulting estimate reads 

(n + n-\m c ug s ta^G^ l5 c\D) « iA<Kg s V% f n f^(0)M 2 D . (23) 

(The amplitude proportional to only the first power of g s reflects the fact that we are not 
yet in the asymptotic heavy quark regime.) 

In the case of photonic 0\ the partonic rate itself describes the probability of a different 
process, D — > 7 + hadrons. Instead we need the similar partonic rate for the photon 
conversion into a ci-quark pair: 

r„=/£ v^)~N.e d , r.,(V) = ^ (l^) 2 ^)^). (24) 

The kinematic integral equals to \n.m 2 c /m 2 — | + the lower cutoff can be taken at 
m around 400 MeV, to match the overall hadronic polarization contribution to charge 
renormalization. Then the integral turns out numerically close to unity. Using this as a 
counterpart to Eq. ( 1221) we get 



(n + Tr-\m c ueia^F^ 5 c\D) ^ i8V2naq d Uf^{0)M 2 D , (25) 

where the small deviation of the explicit log factor from unity has been neglected. 

To cross-check the meaningfulness of the estimates we explored alternative ways. For 
the gluon bilinear O2 this was relating the gluon field operator to the chromomagnetic (or 
kinetic) expectation value in heavy mesons, i.e. treating it as fully nonperturbative. This 
resulted in a different estimate 

2/y 2 

(7r + 7r-\m c ug s ia^G^ l5 c\D) w i-^f^(0)M 2 D . (26) 

Jtt 

Assuming numerically that N c a s ~ 1 we would get a number only 10% smaller than 
Eq. fl23l . The two totally different estimates led to close values because charm lies in 
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between the light- and the heavy-quark regimes where both perturbative partonic and 
nonperturbative description can qualitatively be applied. 

In the case of 0\ photon cannot per se be nonperturbative; instead one can consider 
the photon loop with the gluon emitted internally. This effect is proportional to q u 
and is therefore physically different. It describes the order-a operator mixing of the 
electromagnetic 0\ into chromomagnetic O2. Replacing InAuv by unity we would obtain 

{ti + ti- \m c ueia^F^ lb c\D) « i2y/3 q u ag s Uf^(0)M 2 D . (27) 

Comparing this to Eq. (|25p we find the overall factor different: 

4y/3ira s q u vs. 8y/2irq d , (28) 

with the former 'direct' contribution expectedly dominating since it does not suffer an 
additional perturbative loop factor for gluons. This is partially offset by the larger w-quark 
charge. Altogether the direct photon conversion estimate appears a few times larger, and 
we adopt Eq. fl25|) for the 0\ estimates. 

Collecting all the expressions we arrive at 

, 5.2-icr 2 IT , o.io-icr 3 

IrtlCl ~ 1 — ; ; r , Im C2 



I sin 5 FS i I " |sin5 F si| 

2-10- 3 |T , 4.6-10" 3 

Imc 3 « — — - — r , Imc 4 w — — - — r , (29) 

I sindpsil |sm5 FSI | 

assuming that a particular operator is the sole source of the New-Physics CP violation. 



3 The Neutron EDM 

EDMs in general and specifically the neutron EDM d n are very sensitive probes for physics 
beyond the SM, in particular for CP violation. As already mentioned, a flavor- diagonal 
CP violation with a size of coupling found for the new physics operators in the last section 
would grossly violate the bound, which currently lies at j7] 

\d n \ < 2.9-l(T 26 e-cm. (30) 

In the following we assume that only the AC = ±1 operators induce the non-SM CP 
violation, and estimate their effect. First we recapitulate the salient points of the estimates 
within the Standard Model, where recently a new perspective has been proposed [8]. We 
will discard the possibility of strong CP violation assuming that the long-standing strong 
CP problem will find a solution where the QCD ^-parameter is sufficiently close to zero. 



3.1 The Neutron EDM in the Standard Model 

The estimates of the neutron EDM in the SM have a thirty year long history; the modern 
perspective can be found in review [9]. EDMs may emerge from the second order on in 
the weak interaction and are generally proportional to G\. 
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Motivated by the qualitative success of the constituent quark models in understanding 
of the properties of hadrons, the early estimates of nucleon EDM focused primarily on the 
EDMs of quarks d q . It turned out that for quarks the KM prediction is further suppressed: 
the sum of all the two-loop diagrams vanishes and d q emerge first at the three-loop level 
where an additional loop with at least a gluon must be included [10J. On top of this, the 
quark EDM has to be proportional to the quark mass; this yields an additional suppression 
for the light quarks. The same applies to the color dipole moments of quarks considered 
as the simplest induced CP-odd strong force generated through weak interactions at small 
distances. The unfortunate feature of the quark EDMs is that there is a strong numeric 
cancellation between the leading logarithmic and the subleading terms in the dominating 
EDM of the <i-quark dd [llj . which makes it difficult to to make a definite prediction 
beyond an estimate 

\d u4 \ < 0.5 ■ KT 34 e-cm. 

It has been noted a while ago [121 D2] that the strong suppression intrinsic to u and 
d quarks can be vitiated in composite hadronic systems like nucleons. The transition 
dipole moments changing <i-quark into s-quark, electromagnetic or color, are suppressed 
by the strange quark mass m s , and such flavor- changing transition without a quark charge 
change are already in the loop-induced short- distance renormalization of the bare weak 
interaction due to the so-called Penguin diagrams [I]. 

It is notoriously difficult to account for the long-range part of the strong interactions 
generating the neutron EDM [9]. Usually it is considered that the principal effect comes 
from the diagram in Fig. [TJ having a chiral singularity, or that at least it fairly represents 
the magnitude of d n . One of the vertices in the diagram is a conventional CP-conserving 
AS = 1 weak interaction while the second is the CP-odd Penguin-induced amplitude 
originating from short distances which naturally incorporates heavy quarks, in particular 
top. 




Figure 1: The chirally singular diagram exemplifying the conventional Penguin-based contri- 
bution. One of the vertices is the usual CP-conserving weak amplitude while another contains 
the CP-odd Penguin- mediated operators. 

It has recently been argued [8] that there is a complementary mechanism generating d n 
to the second order in Gp which does not involve short-distance loop effects and is likewise 
free from chiral and SU(3) suppression. It scales like l/m 2 c and would fade out quickly 
for sufficiently heavy charm, yet it may actually dominate d n in the SM since charm is 
marginally heavy in the hadronic mass scale. It originates at the energy scale around m c 
due to interference of the conventional AC = 1, AS = weak amplitudes, much in the 
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same wave as the CP-odd .D-decay asymmetry discussed in Sect. 12.11 Consequently, this 
mechanism would be present, with a modified strength, in the BSM scenarios affecting 
Aacp. The analysis of the BSM contributions to d n presented in Sect. 13.21 parallels the 
SM case, therefore we remind below the main steps of Ref. [8]. 

The observable CP-odd effects appear in the second order in C w and thus are propor- 
tional to Gp, being embodied in 



The generalized GIM-CKM mechanism ensures that the CP-odd piece of £2 is finite in 
the local four-fermion approximation. The conventional form of C w applies to the high 
normalization point around My/. We will neglect, for the most of the consideration, the 
perturbative gluon corrections, since the effect exists even without loops. This makes the 
analysis simpler and more transparent. 

Descending to a low normalization point we first integrate out top quark and at the 
second stage, below also the bottom quark. At tree level integrating them out simply 
means discarding all the terms containing t- or 6-quark fields. As a result, we arrive at 
£2 generated by a superficially two-family weak Lagrangian 



In fact, this is not a true two-family case, since the four Vu do not form a unitary 
matrix; in particular, it is not CP-invariant. The phases in the four remaining CKM 
couplings cannot all be removed simultaneously by a redefinition of the four quark fields, 
as quantified by A in Eq. (Q. 

Interested in d n we need to consider only the terms in £ 2 that are diagonal in all 
four quark flavors. Moreover, we can omit explicitly CP-invariant terms of the form 
of a product of an operator with its conjugated. Only two operators remain after this 
selection, those proportional to the CKM product V* s V us V c dV* d , and also their Hermitian 
conjugated. These non-local 8-quark operators include both q and q fields for each of the 
four quark flavors: the CP-odd invariant A (as well as CP-violation altogether) vanishes 
wherever any single CKM matrix element becomes zero. The two terms in £2 differ by 
the type, up- or down-, of the quark-antiquark pairs coming off the same weak vertex, see 
Fig. [2J The CP-odd amplitudes conventionally considered for d n are of type a) where one 
of the weak vertices has cc and another uu, and both have AS' = —AD = 1. The type b) 
amplitude has the dd and ss pairs in the two weak vertices, respectively, while each has 
|AC| = 1. These were routinely omitted. 

Considering the nucleon amplitudes we need to eventually integrate out the charm 
quark field as well. At this point the distinction between the two types of terms becomes 
important. Where the two charm fields belong to the same four-fermion vertex in the 
product Eq. (}3~Tj) . Fig. [2^, they can be contracted into the short-distance loop yielding, 
for instance, the usual perturbative Penguins. These are the conventional source of the 
long-distance CP-odd effects [121 H3]- The loop cannot be formed for the alternative 




(31) 



■w 



Aj^ with J ' = V cs cT.s + V cd cT„d + V us uTuS + V ud uTud . 



(32) 
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possibility where c and c belong to different £ w , Fig. [2b since the charm quark must 
propagate between the two vertices. This is the reason why such contributions were 
usually discarded. 




Figure 2: Two types of CP-odd terms. Weak vertices must be off-diagonal in flavor, either 
for down-type (a) or up-type (b) quark. Solid dots denote the four-quark vertices. Light lines 
correspond to u, d or s quarks, thicker lines stand for charm. 

Nevertheless, the latter term has an advantage: it does not involve short- distance 
loops, and has a single charm propagator, although highly virtual in the hadronic scale. 
Each weak vertex contains a flavorless quark-antiquark pair, but these are light down-type 
quarks d and s and are not contracted via a perturbative loop; they will go instead into 
the nucleon wavefunction. The corresponding operator is 

V cs V; d V ud V: s J d A x zT{(dT>)(?TV)o ■ (cT v s)(sr v u) x + H.c. (33) 

The Hermitian conjugate, apart from complex conjugation of the CKM product, is simply 
the exchange between s and d, s^-d (this particular property does not hold beyond the 
SM). For the sake of transparency, we have passed here to the sum and the difference of 
the two operators in C w , in terms of Eq. (@|. 

As the space separation x in Eq. f )33|) is of order 1 /m c , eliminating charm results in a 
local OPE; the expansion parameter /ihadr/^ic is not too small and we need to keep a few 
first terms. The tree-level OPE is particularly simple here and amounts to the series 

c(0)c(x) = ( 1 —) = —5\x) + -^5\x) tip + \d\x) {tip) 2 + ... (34) 

\m c -i]p J 0x m c m z c m% 

valid under the T-product. For purely left-handed weak currents in the SM the odd 
powers of l/m c in Eq. ( 1341) are projected out, including the leading l/m c piece. We then 
retain only the term and arrive at the local effective CP-odd LagrangiarH 

£_= -tA-^d uds , (35) 

O uds = {uT^d) (Jr M 0r^) {sT»u) -{d^s} = 

(uT^d)-[(dT^iIpT„s)-(sT u u) + (dT^ a T u s)id a (sT u u)] - {d s}; 

2 We have changed notations compared to Ref. [8]: now O u ds, O u ds and 0" ds all include subtraction 
of the Hermitian conjugated operator. 
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in the last expression the covariant derivative acts only on the s-quark field immediately 
following it. 

To address the electric dipole moments we need to incorporate the electromagnetic 
interaction. One photon source lies in the covariant derivative in the operator O u d s , which 
includes electromagnetic potential along with the gluon gauge field. It is proportional to 
the up-type quark electric charge +|. The corresponding photon vertex is local and is 
given by the Lorentz- vector six-quark operator which we denote as 0® ds : 

0: ds = (u^(l- l5 )s) (s 7 ^7 a 7,(l-75)d) (d7"(l-75)w) - {d s}. (36) 

Another, non-local contribution is the T-product of the pure QCD part of O u d s 

O uds = (uY(l-l 5 )s) {s ltl lv {l- l5 )d) (d 7 "(l-7 5 )t0 - {d s} (37) 

with the light-quark electromagnetic currently The total photon vertex is thus given by 
the effective CP-odd Lagrangian 

C = El?^ ( 38 ) 
q 

where is the electromagnetic potential and e is the unit charge. 

In principle, the local and non-local pieces above correspond to distinct physics: one 
has photon emitted from distances of order l/m c while the latter senses long-distance 
charge distribution. The latter usually dominates, however the specifics of the left-handed 
weak interactions in the SM make them of the same order. 

An important feature of the considered contribution is that it remains finite in the 
chiral limit and it does not vanish if d and s quarks become nearly degenerate, at first 
glance contradicting the origin of the KM mechanism where an additional SU (2) freedom 
to mix s and d makes the theory CP-invariant at m s = rridQ This in fact is fully consistent, 
since the external state, the neutron, is explicitly s -H- d non-symmetric, and would not 
stay invariant under the mixing transformation. The same applies, for instance, to the 
(i-quark EDM. In contrast, short-distance effects of light quarks in the loops involve 
severe GIM-type suppression proportional to the powers of the light quark masses. The 
quark electric dipole moments as the purely short-distance contributions are explicitly 
proportional to the corresponding quark mass for chirality reasons. 

Let us parenthetically note that there is no a formal contradiction between the non- 
vanishing expression for d n and the explicit T-invariance at literally m s = rrid either. 
T-invariance prohibits dipole moments only for the eigenstates of the Hamiltonian; for 
instance, non-diagonal dipole moment matrix elements are perfectly allowed. The above 

3 In general only the sum of the two terms yields the transverse electromagnetic vertex; however, when 
projected on the dipole moment Lorentz structures they separately conserve current. 

4 The contribution would vanish if charm and top become degenerate; considering the cases of de- 
generate bottom and strange quarks, or charm and up makes no sense in this context since it has been 
assumed as the starting point that mb,m c 3> /ihadr while u, d and s are light quarks. 



A n C„ 



G 

mi 



o: ds + d*xtT{o uds (o) j: m (x)} 
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considered neutron, a baryon state with strangeness S = is such a physical eigenstate 
only as long as the mass splitting m s — remains much larger than the weak corrections 
to the hadron masses, since weak interactions violate flavor. Where m s — md becomes of 
the order of GFm q w? c mixing between neutron and its strange partners must be accounted 
for. In this regime the time-violating EDM for a physical state should be distinguished 
from the conventional d n ; the former would vanish at m s = rrid- 

The CP-odd operators contain strange quark fields. This means that the induced 
effects would vanish in a valence approximation to nucleon where only d and u quarks 
are active. It is known, however, that even at low normalization point the strange sea 
in nucleon is only moderately suppressed. The large-iV c perspective on the nucleons 
paralleling the picture of the baryon as a quantized soliton of the pseudogoldstone meson 
field [2] makes this explicit: the weight of the operators with strange quarks in the chiral 
limit is generally determined simply by the operator-specific Clebsh-Gordan coefficients 
of the SU(3) group. 

Such an 'intrinsic strangeness' suppression is specific for the considered mechanism 
to generate d n in the SM; the conventional contribution trades it in for the 'intrinsic 
charm'. Associating the virtual-pair suppression with the strangeness sea in the nucleon 
is probably a relatively light price to pay. In contrast, the perturbative Penguin effects 
yield small coefficients whenever considered in the truly short-distance regime. 

3.1.1 Matrix elements 

The CP-odd operators 0^ ds and O u d s have high dimension; this is routinely associated 
with being poorly defined for practical applications. However, these particular operators 
possess intrinsic symmetry properties, including antisymmetry in respect to s -<->■ d, which 
prohibit mixing with lower- dimension operators, and make them a suitable object for the 
full-fledged nonperturbative analysis. 

The neutron EDM is obtained by evaluating the hadronic operators in Eq. fl38|) over 
the neutron state. Since is T-odd, the matrix element vanishes for zero momentum 
transfer and the linear in q term describes d n : 

(n(p+g)|£^|n(p)} = d n q v u n {p+q)ia> a '75U n {p). (39) 

Neither of the two matrix elements involved are easy to evaluate, although one may 
hope that such a contribution may eventually be determined without major ambiguity, 
including the definitive prediction for the overall sign. Although only the P-violating part 
of O^J s contributes, the original form is more compact and makes symmetry explicit. 

The contact operator 0^ ds is a product of three left-handed flavor currents; O u d s 
instead of the sd current has a flavor non-diagonal left-handed partner of the quark energy- 
momentum tensor in the chiral limit. Therefore it seems plausible that the required matrix 
elements can be directly calculated within the frameworks like the Skyrme model [151 EH] > 
or in its dynamic QCD counterpart [16] derived in the large- N c limit from the instanton 
liquid approximation. 
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Lacking presently more substantiated calculations we resort to the simple dimensional 
estimates. For the local piece we put 

(n(p+q)\OZ ds \n(p)} = 2iJC uds q v u(p+q)i<r ,a ''y B u(p). (40) 

The reduced matrix element )C u d s has dimension of mass to the fifth power. We estimate 
it as 

I^QwisI ~ K A'-hadr) (41) 

where /ihadr is a typical hadronic momentum scale and k stands for the 'strangeness 
suppression' to account for the fact that neutron has no valence strange quarks; re«l/3 
is taken typical guess. 

Due to the high dimension of the operators the estimate for d n depends dramatically 
on the value used for /ihadr- Although the typical momentum of quarks in nucleon is 
around 600 MeV or higher, using this as /ihadr would strongly overestimate the effect. Six 
powers of mass in Eq. (14~TI) come from the product of two local light quark currents each 
intrinsically containing factors N c /8ir 2 when converted into the conventional momentum 
representation. This is illustrated by the magnitude of the vacuum quark condensate 
where such a factor effectively reduces //j^dr down to ~ (250 MeV) 3 . 

To account for such differences we assign a factor of (0.25 GeV) 3 = ^ to each additional 
quark current in the product, while the remaining dimension will be made of the powers 
of /ihadr taken around 500 MeV. Then this contribution to d n becomes 




where A~3.4-10~ 5 has been used. An independent enhancement factor may come from 
summation over the Lorentz indices in the currents, but we neglect it. 

The most naive estimates for the d n induced by the non-local piece in Eq. (!38|) would 
yield a similar dimensional scaling except that no explicit factor e c = | appears: the 
dimension of the non-local T-product is the same as of 0^ ds itself. Following the more 
careful way advocated above where we distinguish the mass scale associated with the local 
product of the quark fields, the result is literally different: 

« e a|32 K ,>- » 1.2 • 10-e.cm x K (^_)' (^X); ,43) 

numerically the difference is not radical, however. 

Alternatively, the non-local contributions can be analyzed focusing on the contribu- 
tions of the individual intermediate states, usually the lowest in mass. Among them the 
hidden-strangeness states, including KA(E) look promising suggesting a way to dynam- 
ically estimate the 'intrinsic strangeness' factor k. In the standard model, however, the 
corresponding loops are not infrared-enhanced and rather saturate at large virtual mass 
yielding a result strongly dependent on the assumed cutoff. For the same reason the kaon 
and the lowest baryon as the intermediate state are not any more remarkable a priori 
than ordinary resonances. 
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Ref. jB] considered the contribution of the lowest resonant state, the | nucleon reso- 
nance iV(1535) referred to below as N, as an alternative estimate of the non-local piece 
in d n . In terms of the two hadronic vertices, 

(n{p')\J^)\N(p)) = -PNU n ta^ l5 q uUff , (JV(p , )|O wfe (0)|n(p)) = 16iAU%M n (44) 
the sum of the Feynman diagrams in Fig. [3] gives 

The electromagnetic vertex estimated from the measured transition iV — >n + / y becomes 
(0.34 ± 0.08) GeV -1 . The induced weak CP-odd vertex is estimated in the dimen- 
sional way, for the dimension-ten operator U( i s yielding 

Wuds\ ~ K /ij/ihadr- (46) 

Finally this estimate reads 

w*- A f ^• wr31 -" (^)'yiv) ■ <«» 

This value is consistent with the direct dimensional estimate of the non-local contribution, 
in particular considering the fact that the lowest excited state alone may not necessary 
saturate it. Therefore, in further applications we generally follow the more straightforward 
estimates paralleling Eq. ( 143]) . 



N 



N 



Figure 3: Nonlocal contribution to d n with the intermediate N. Solid block denotes the 
CP-odd operator O u d s - 

Finally, our estimate for d n in the SM centers around 10 _31 e-cm although even the 
values 5 to 10 times larger may not be excluded. 

The natural benchmark for the CKM d n in the SM evidently lies about 

4ocAG^ adr (48) 
In the same terms the loop-less contribution considered above is 

d n oc A G^/ihadr " ' K - (49) 
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The last factor reflects the absence of valence strange quarks in the nucleon. The related 
suppression is unavoidable for d n in one form or another; it is natural to think that 
paying the price by the soft strangeness content in the nucleon state is the minimal burden. 
Therefore, from this perspective such a contribution appears to bear a mild model-specific 
suppression, since m c in practice only moderately exceeds the characteristic hadronic scale 

At the same time d n does not contain parametric chiral enhancement, ln/i^ adr /m^ or 
numerically significant scalar matrix elements possible in the case of generic couplings. 
This implies also a loss of a potential factor of a few; it may be recovered by the interactions 
generated beyond the SM. 

Throughout the long history of the conventional long-distance contributions to d n in 
the SM it has usually been considered [9] that the principal effect comes from the diagram 
in Fig. [1] peculiar by showing a chiral singularity. It has been calculated in Ref. [13], 

rf „ KeG |A%^ln42l<|^i(2 a -l) 9 ,ln^ (50) 
27V27T 3 m A c f n m s m n 

(A is a strong constant parameterizing the conventional CP-even vertex and a a dimen- 
sionless ratio of two £77(3) meson-to-baryon axial couplings, while C pe rt stands for addi- 
tional perturbative factors). The original authors' estimate was close to <i n ~2-10 _32 e-cm. 
It was done in 1981 when even the size of the CKM admixture of the third generation was 
unknown and was thought to be of the scale of 9c- The equivalent of the CP- violating 
parameter A likewise was estimated assuming m t « 30 GeV, yet a value only 1.5 times 
larger than known today, see Eq. ([6]), was used. At the same time the used log ratio of the 
t and c quark masses was somewhat smaller. The size of this contribution is now usually 
cited as d n ~10~ 32 e-CHi [Sj. 

The d n value Eq. (1301) is proportional to a s /Ti from the short-distance Penguin loop. It 
also contains a factor m^ocm 5 compared to the benchmark Eq. (jUJ), however the overall 
light-quark mass scaling is remarkable: m Uj d, enter divided by m s rather than by /^hadr- 
Therefore, in the SU(3) chiral limit where all m Ut d tS — > 0, m q /m s fixed it would stay finite. 

In practice, however the SU (2) chiral limit m u ^ — > 0, m s fixed is more relevant in 
numeric estimates. In this case Eq. ( 150]) has an additional factor of the light-quark mass 
m q compared to the benchmark value Eq. (148]) . This is in agreement with the general fact 
stated in Appendix |U the contribution contains a chiral log and therefore must include 
an explicit factor of m q since the SM weak amplitudes do not contain right-handed light 
quark fields. As emphasized there, it is sufficient to check this for the bare weak vertices. 

This illustrates the underlying problem in estimating d n in the SM: the physically 
distinct chirally singular contributions have to be m g -suppressed. The leading-m g con- 
tributions are not related to soft pions and are rather saturated at the loop momenta of 
the typical hadronic mass scale /ihadr, or by resonances with a significant mass gap. Such 
effects are generally uncertain and may involve cancellations. 

The conventional SM contribution Eq. ( 150]) . therefore, has an additional light-mass 
suppression cx m q on top of the perturbative short-distance factor. Although it is partially 
offset by numerically large factors accompanying the amplitudes with right-handed light 
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quarks, together with the perturbative loop factor it results in a certain suppression. This 
may explain the larger number for the loopless EDM which we estimate to be around 
10~ 31 e-cm. 



3.2 Neutron EDM and a BSM Charm CP Violation 

In order to estimate the effect of the |AC| = 1 amplitudes on the neutron EDM we should 
replace one of the two C w in the product £2 by the New Physics operators. We will assume 
that we get a reasonable estimate when considering one operator at a time; this basically 
corresponds to the assumption that in the neutron EDM we do not have a destructive 
interference absent from the charm decays. 

As is clear from the analysis of the SM case, the operator structure obtained upon 
integrating out charm depends on its chirality in the NP amplitude: in the left-handed 
case it follows the SM case. Where the charm field is right-handed, only the odd-power 
terms l/m c , 1/m 3 ... survive. In this case the leading term suffers less from the fih&dr/ m c 
suppression, however it does not include the leading contact photon vertex (the photon 
operator 0% is an exception in this respect) which appeared to yield a few times larger 
contribution, at least within our estimates. The contact photon vertex is then delayed till 
order l/m? c . Such a peculiarity introduces certain difference, but in view of the relatively 
mild numeric power suppression the presence of the nonlocal T-product term to the leading 
l/m c order for the right-handed charm does not appear to bring in a notable numeric 
difference. 

The chiral content of the light valence quarks generally makes a bigger difference. 
The nucleon matrix elements with both left-handed and right-handed fields are usually 
numerically enhanced as seen on the example of the nucleon cr-term. Moreover, a CP-odd 
scalar pion-to-nucleon coupling may be induced. Although the CP-nonconservation case 
is more involved, it can be stated that this vertex at small momentum transfer would be 
proportional to the light quark masses and therefore negligible in practice unless the New 
Physics operators include right-handed light quarks, see Appendix HI 

The contact photon vertex contribution to d n does not depend on the induced pion- 
nucleon interaction. The scalar pion vertex Q s pmr~, on the other hand, generates a 
chirally enhanced long-distance contribution to the T-product piece with the ir~p inter- 
mediate state, described by the diagrams in Fig. HJ 

GsQnNN , A 2 _ . Qs9a , A 2 _ „. , K1 x 

with A the ultraviolet cutoff. In the actual world the chiral ln constitutes a moderate 
factor about 3 and therefore is remarkable more in the conceptual aspect. However it 
comes proportional to the large couplings and this makes up for the loop factor. In 
the few considered examples the overall chiral enhancement roughly offsets the typical 
suppression (at the same order in l/m c ) of the T-product piece relative to the contact 
photon vertex contribution. Regardless of the details, it is clear that the chiral log per se 
is a too weak singularity to change dramatically the expected magnitude of d n . 
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Figure 4: The chirally singular diagram generated by the scalar vertex. One of the vertices 
is the usual strong CP-conserving pseudoscalar coupling while another is the induced CP-odd 
scalar vertex. 



The effect of the scalar pion-to-nucleon vertex can be more pronounced in atomic 
EDMs. The pion-mediated nuclear forces are relatively long-range and may be addition- 
ally amplified for the isoscalar coupling in heavy atoms. 

The current algebra technique allows to unambiguously determine the induced scalar 
vertex in the chiral limit, see Appendix HI Of course, in the general case it would require 
the nucleon expectation values of similar high-dimension effective operators, currently 
estimated in a rather crude way. An additional uncertainty would come from possible 
cancellations due to the second, pole-subtraction term in Eq. (lA.lj) . As a rule, aiming 
only at the overall magnitude, we simply neglect this term. 

Contracting charm propagator we end up with the multi-quark CP-odd operators. 
The neutron EDM induced by them is evaluated applying the dimensional estimates 
elaborated for the SM case; they are described in Sect. 13.11 The non- valence strange 
quarks are neglected here and no factor k appears. 

Considering the quark bilinears, we start with the more natural gluonic 2 . In view 
of the above mentioned difference in the OPE, we consider separately the cases of 
and O2 containing right- and left-handed c-field, respectively, rather than its scalar and 
pseudoscalar versions. 

We start with 0<f = m c cg s iaG{l — ^)u. Here charm induces already the l/m c effect 
generating d n via the T-product with the electromagnetic current. The contact term 
emerges only to order 1/m 3 , and we discard it. The corresponding CP-odd operator is 

O(-) = -zImc 2 ^ Sin2g cCQs2gc [ug s ia a eG a ^(l- l5 )ddj»(l- l5 )u- s] + H.c. 

(52) 

For the scale estimate we simply factor out the operator g s ia a pG a ^ and assume it has the 
value similar to the one in heavy mesons or baryons, 2/j,q 0.7 GeV 2 . A close value is 
obtained if we use the vacuum condensate (qg s iaGq) ~ 0.8 GeV 2 (qq) . Discarding strange 
quarks and applying to the rest our dimensional estimate we get 

(0 25 GeV) 6 

\d n \ « Imc 2 G! sin 2 # c cos 2 fl c 32/4^- rip- Xfl ~ 1-1 ■ 10" 26 e-cm ■ Im c 2 Xfl , (53) 



(/^hadr 



where xa ~ 1 to 2 is a flavor factor which accounts for the fact that there are two d quarks 
in the neutron. Using Eqs. (1291) for Imc 2 we get 

\d n \ « 10~ 30 , , ^ fl re-cm « 2.3xa • 10~ 30 e-cm . (54) 

SHI Oir+ir- 
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Here and in what follows we assume | sin5psi| ~ 0.5 as a typical value. Thus we expect 
an enhancement of roughly a factor of thirty. 

Now we turn to 0\ —m c c.g s iaG{l-\-^^)u. Here the leading term l/m c vanishes like in 
the SM and the l/m c expansion starts with 1/m 2 , yet we have the right-handed w-quark 
which entails chiral enhancements in the nucleon matrix elements. There are two distinct 
contributions to the leading order like in the SM, the contact and non-local. For the latter 
we obtain 

H 2G^sin 2 g c cos 2 fl c 

( '= -%\mc 2 — : x (55) 

m\ 

% s fc^-i(^+<^V a „-^ a ^ - d ^ s] +H.c. 

and the contact electromagnetic vertex is given by 
n»(-) -t 2G 2 ,sin 2 # c cos 2 fl c 

( '= -dmc 2 — : x (56) 



m c 



% s (G^-i(>"+G^ a „-^ a G^ -d^s 



-H.c. 



It likewise may have a numeric chiral enhancement due to right-handed w-quark, yet no 
literal chiral log from the pion loop. Our estimate for it reads as 

\d n \ loc « ^Im c 2 G% sin 2 # c cos 2 # c (0-25 GeV) ^ ^ ~2.5 • 1(T 26 e -cm • Im c 2 x& Xscsi, 

(57) 

where we, as above, have equated the whole bracket containing the gluon field strength, 
including g s , with 2/^; yet another factor x sca i has been added to indicate a possible 
enhancement of the scalar expectation value (cf. the size of the nucleon cx-term). This is 
about 3 times larger than in Eq. f l5"3"|) . 

The non-local contribution estimated dimensionally is typically 2.5 to 3 times smaller 
than the contact one. However here the right-handed w-quark induces the nonvanishing 
scalar pion-nucleon vertex and the 7r~p intermediate state yields a chiral log, cf. Eq. (!5T|) . 
Combined with the current algebra result for the scalar version as described above this 
enhancement numerically turns out about 3.5, i.e. we get a number close to the contact 
estimate Eq. ( |5"T|) . 

Thus, we can use for this case the local estimate Eq. (1571) and Imc2 from Eq. fl29l) . 
Then 

|41 sa 2.5-HT 30 . Xfl f cal e-cm « X R Xscai ■ 5- l(r 30 e-cm . (58) 

| Sill Ott+tt- I 

For this chiral structure we get about an 80-fold enhancement compared to the SM. 

The photonic operators 0\ are the simplest case since only the contact photon vertex 
should be considered to the leading order in a. For the case of the right-handed c quark, 
0^ = em c ciaF(l— , ~f 5 )u the leading term in the l/m c expansion yields 

0" (-) = 2zImciG 2 7 sin 2 fl c cos 2 fl c d M [uia fiv j a (l-j5)dd^ a (l-j 5 )u - d s] + H.c. (59) 
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The matrix elements of the CP-even partner of such an operator may have been estimated 
in the literature. Applying our standard recipe we get 

|d n | « 16|Imci|G|sin 2 ^cos 2 ^ c (0.25GeV) 3 Xfl 

« 3.2 ■ 10" 26 e-cm- ImcrXfl ~ 3.4-10~ 27 e-cm-Xfi . (60) 

Owing to its nature this operator yields the largest enhancement of all the new physics 
operators. Nevertheless it is still safe in respect to experimental bounds. 

The operator 0\ with the opposite chiralities, Of = em c caF ( 1 + 75) u has a mild 
suppression by a factor ^thadr/^c, however it can be enhanced by larger matrix elements 
appearing with the right-handed u quark. Therefore we expect to have here the same 
numeric estimate as for Of, within a factor of 0.5 to 2. 

Finally we consider the four-quark operator O4. This is the case of both the leading- 
order 1 / m c contribution and of the chiral enhancement from the light valence quarks in 
the nucleon. Neglecting the strange quarks we have 



O(-) = nmc 4 ^^ [^(1-75)^^7,(1-75)^^(1-75)^] +H.c. (61) 

m, 



The contact photon interaction would come suppressed by two powers of l/m c . On the 
other hand, the above leading-m c interaction enjoys a chiral pion loop enhancement in 
the T-product with J cm . Taking the axial charge commutator and neglecting the pole 
subtraction term in Eq. (lA.lj) the scalar pion vertex becomes 

G 2 sin 2 cos 2 

Q s UpU n = -ilmc 4 — 2{p\u^(l-j 5 )dd^^ u (l+7 5 )ddY (l-l5)d\n) g =Q . (62) 

According to our dimensional rules this amounts to 

\r \ it 1 fi G>sin 2 c cos 2 c , , 6 2 , . 

\g s \ = Imc4 8 (0.25GeV) D XscaiX fl (63) 

m c f n 

and results in 

Kl « |Imc 4 | ° l S[n Z C COS26c 9 a In % (0.25 GeV) 6 Xs cai xl 

^ 6-10 _28 e-cm- Imc 4 XscaiXfl- (64) 
Using the estimate Eq. ( )29|) we end up with 

|4| ~5.7-10- 30 e-cm Xs caiXfl- (65) 

Should we apply the simple-minded dimensional estimate to the T-product contribution 
without considering specifically the pion loop or paying attention to the potential chiral 
enhancement, we would get a somewhat smaller but a consistent value 

\d n \ w 4-10" 28 e-cm- |Imc 4 | xl ~ 4-10" 30 e-cmx 2 1 . (66) 
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Table 1: The estimated D° decay amplitudes, the strength of the CP-odd couplings and ex- 
pected d n . The two sub-columns for the chromomagnetic operator O2 correspond to the left- 
handed (left) and right-handed (right) charm fields, respectively. 
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Therefore, in the case of O4 the induced d n is about 100 times the SM. The origin 
is evident: O4 has a color structure unfavorable to D — > 7t + tt~ . At the same time, the 
chirality choice is optimal for both charm and light quarks, in the nucleon matrix elements. 
The combination of the two yields an additional factor of 10 enhancement in our estimates. 

For convenience, Table [T] summarizes our estimates of d n in this section along with 
the values of Imc^ from Sect. 12.21 

3.3 A comment on the atomic EDMs 

The atomic size exceeds the nucleon radius by several orders of magnitude, and as a 
matter of principle they may have larger EDMs; in particular, this applies to paramagnetic 
atoms. The enhanced EDM, however may originate there mainly through T-violation in 
the lepton sector, with the electron EDM itself or via the induced contact interaction with 
the nucleons. Such manifestations of New Physics are not directly associated with the 
milliweak interaction of quarks and are beyond the subject of the present study. 

In diamagnetic atoms like mercury the screening mandated by the Schiff theorem is 
rather effective and the overall EDM appears to be dominated by the induced isoscalar 
CP-odd ir°NN coupling affecting non-pointlike electromagnetic potential of the nucleus 
- yet still at a rather suppressed level, 

d Hg « £ s ■ 3.5 TO" 18 e-cm, (67) 

see Refs. [T71 [9]. Using, for instance the estimate Eq. fl63l we can expect for the isoscalar 
coupling \Q S \ ^10~ 15 . Therefore, as anticipated the diamagnetic atom EDMs, while prob- 
ably not yet fully competitive in sensitivity with the direct d n , may become comparable 
in certain NP scenarios yielding amplitudes with a right-handed light quark, owing to the 
recent radical improvement [18] in the precision for the 199 Hg EDM. 
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4 Conclusions 



The KM mechanism of CP violation in the Standard Model is an instructive example of a 
realistic phenomenological theory where the dominant contribution to the electric dipole 
moment of neutron comes not from the effective CP-odd operators of lowest dimension, 
but via a nontrivial interplay of different amplitudes at a relatively low energy scale. In 
the SM this evidently roots in an intricate nature of the CP violation intimately related 
to flavor dynamics requiring existence of a few generations. It may be interesting to 
investigate, in general terms, if a similar pattern can naturally fit theories beyond the SM 
which would describe flavor dynamics at a more fundamental level. 

We have argued that in the SM itself with vanishing #-term the neutron electric dipole 
moment has natural size about 10 _31 e-cm and may even exceed this, due to the interfer- 
ence, at the momentum scale around 1 GeV of the two AC = 1, AS' = weak four-quark 
amplitudes. This mechanism does not require short-distance loop effects, is finite in the 
chiral limit and does not depend on the strange- vs. down-quark mass splitting. 

The CP-odd direct-type D° decay asymmetry reported recently at the level of 10~ 2 
does not naturally conform the expectations in the SM, which are typically an order of 
magnitude smaller. This may be an indication for new CP-odd forces beyond the SM, 
although such an interpretation should still be viewed cautiously. 

If New Physics indeed induce a milliweak CP-odd decay amplitude in charmed par- 
ticles, it may also be expected to generate, at the NP scale, flavor-diagonal CP-odd 
interactions in the light hadron sector. The EDMs of nucleons and atoms are extremely 
sensitive to them, and the existing experimental bounds place strong constraints on the 
effective interactions seen at low energies. Such low-energy effective operators are model- 
dependent and their connection to the charm CP violation is indirect, to say the least. 

Nevertheless, a certain, possibly subdominant contribution to d n is generated at the 
charm energy scale in a direct analogy with the Standard Model. It is fully independent 
of the effects originating from the NP scale and directly reflects the scale of CP violation 
in charmed particles. Our analysis suggests that this would increase d n compared to 
the SM prediction by more than an order of magnitude: the typical enhancement is 
between 30 and 100, depending on the chiral, color and flavor composition of the charm 
NP amplitudes. In an ad hoc case of the CP violation through the electromagnetic c — >u 
dipole operator alone the neutron EDM can be even as larger as 5-10~ 27 e-cm. However, 
the possibility itself for NP to generate such a CP-odd electromagnetic operator but not 
a similar chromomagnetic one of the commensurate strength, does not look natural. 

We conclude that New Physics CP violation in charm at the reported level remains 
safe in respect to existing strong experimental bounds on EDMs, as long as the direct 
effects are considered. At the same time it would significantly reduce the gap between 
the bounds and the expected size of the EDMs, and would make the new generation of 
the EDM experiments more topical. 

In the present analysis we have assumed that a new source of CP violation appears 
solely in |AC| = 1 interactions. Eventually the known flavor dynamics must be embedded 
in a full picture of flavor together with CP violation at some high scale, where new 
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dynamic fields are also present. Attempts to investigate the new phenomena along these 
lines have been reported in [5JE], considering the observed Aacp e.g. in a supersymmetric 
framework. This generically induces additional CP violation compared to our scenario, 
which would modify the impact on the EDMs. 
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Appendix: Generalized Goldberger-Treiman relation 
and the scalar ttNN vertex 

CP-odd perturbations in general induce the parity- violating scalar pion-to-nucleon vertex 
which may not vanish at small pion momentum. Such amplitudes often play a special role 
owing to the small (vanishing in the chiral limit) pion mass. Here we give a compact cur- 
rent algebra derivation of the corresponding small-momentum limit for a general CP-odd 
perturbation operator 0~ . We also point out that the induced CP-odd vertex necessar- 
ily vanishes in the chiral limit at zero pion momentum for any operator 0~ which does 
not involve the right-handed light quarks (their absence may be established at arbitrary 
chosen normalization point). 

For simplicity of the notations we consider the charge pion. Its amplitude off the 
nucleon takes the following form in the limit of vanishing pion momentum: 



(N\imo),oMUN) - <0|llQg '° ) ;"l;' lo|0> ww»\d 

-(0\il)ip\0) 



(A.l) 



where Q§ — u^^d, = ^y-^ 5 '? are ^ ne ax ^ charge densities and the equal-time commu- 
tators marked with the null subscript are calculated according to the standard rules. 

To prove it, we first establish a counterpart of the Goldberger-Treiman relation for 
the general case where parity can be violated. To this end we consider the exact nu- 
cleon matrix element of the non-singlet light-flavor axial current J M (let it be Jj" 5 , for 
concreteness) 

(N\J*\N) = g A (q 2 )^Nl^N + b(q 2 )^> N a^q u ^N + C(q 2 )q^ N ^ N + 

a(q 2 )^! Nl ^ N + b(q 2 )^ N a^q v ^ N + c(q 2 )q^ N ^ N , (A.2) 

where the last three terms violate parity being induced by Consequently, 

(N\d^\N) = (2M N g A (q 2 ) + q 2 C (q 2 ))^ N ij^ N + ic{q 2 )q 2 ^ N ^> N . (A.3) 
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which gives two relations, for the pseudoscalar and for the scalar structures. 

Noether theorem relates the divergence of the current obtained from the quark equa- 
tions of motion to the variation of the Lagrangian under the chiral symmetry transforma- 
tion; in the case of QCD the variation comes from the conventional light quark mass term 
and from the corresponding commutator of the axial charge with which we denote 
by -W+E 

d,Jl = 2m q m 75 d + V + , = f[Q+(0),O(-)] , (A.4) 

and therefore 

(N\2m q m l5 d + T>M\N) = (2M N g A (q 2 ) + q 2 C(q 2 ))^ N i lb ^> N + ic(q 2 )q 2 ^ N ^ N . (A.5) 

This equation can be taken in the limit m q — t-0 and to the first order in the 0^> pertur- 
bation. The pseudoscalar structure becomes the Goldberger-Treiman relation f^g^NN — 
2M N g A stating the existence of the Goldstine boson through the pole in C(q 2 ) as long 
as M N g A {0) 7^0. The scalar term dictates that the pion pole residue likewise carries the 
scalar nucleon vertex proportional to {N\D^ \N) at zero momentum transfer: 

(n\vM(o)\ny 

A tNN = . (A. 6) 

JTT 



X>( + ) in Eq. (1A.4|) is just the conventional PCAC commutator. The subtlety is important, 



however that the matrix element above stands for the exact nucleon states rather than 
for the unperturbed QCD ones as is usually implied when expanding in perturbation; 
this fact is indicated by the prime in Eq. flA.6p . The difference becomes important in the 
chiral limit where the pion mass is parametrically small, as illustrated later. 

To bypass this complication we apply a Lagrange multiplier trick, namely consider, 
instead, the CP-odd perturbation ^ =0^ — \(ui~f 5 u—di'y 5 d) with an arbitrary A, and 
keep m q nonzero. On one hand, the operator ui^u — di^ 5 d = is the total 

derivative in QCD and does not change any strong amplitude; hence it can be added to 
the perturbation for free. On the other hand, A can be taken such that the perturbation 
O x becomes nonsingular in the chiral limit. (In other words, in this case one can safely 
perform a double expansion in m q and in 0^~\) Since the chiral singularity comes from 
the pion pole in the correlators, the value of A is determined by vanishing of the residue 
WO'-'lO}, or 

a = . im ,, <°i°'-v> = lim mmw,Q'-^ 

™ q ^o (0\ui-f 5 u-di^d\TT°) ™ 9 ->o -2(0|^|0) 

With this choice of 0^ x ^ the exact nucleon states in Eq. flA.61) enjoy a regular expansion 
in both m q and in free from a l/m q enhancement. Therefore, to the first order in O x 
at m q /ihadr w e can ignore the difference between the exact and the unperturbed nucleon 
states in Eq. (IA.6I) . The net effect of the resummation of the pole terms then amounts to 



3 For simplicity we assume m u —md- 
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subtracting from X>( + ) the divergence of the current with the mass term X^uij^u — dij^d) 
which is equal to 2Xud. This is the relation Eq. (lA.ip . 

The presence of a subtraction term beyond the conventional PCAC commutator in the 
weak transition amplitudes has been appreciated in the context of electroweak calculations 
in the early 1980s [191 120], being evident when using a- models to visualize the chiral 
symmetry breaking. As noted in the end of this section, it is likewise intuitive in the 
ordinary CP-conserving weak decays, in particular of kaons, without recourse to chiral 
Lagrangians. Later it was more systematically incorporated in the chiral expansion for 
many EDM calculations beyond the SM [2T] . 

The above proof, while short and general, may look somewhat mysterious since a finite 
yet calculable part of the usual PCAC commutator term appears to be miraculously eaten 
up only as a result of the failure of the conventional chiral expansion. The cover of mystery 
is removed once the corresponding diagrams are identified and are accurately calculated. 
This is possible using the double expansion, in and then in m q . We illustrate this in 
what follows. 

The pion-nucleon amplitude to the first order in perturbation has two pieces given 
by the irreducible and the pole diagrams, respectively, see Fig. [5j The latter are those 
which become singular in the limit of m q — > or at vanishing pion momenta. We need to 
consider them in the kinematics where n~ has a finite momentum yet small compared to 
the hadronic scale /ihadr, while 7r° has nearly vanishing momentum driven down by small 
m q . The pole diagrams in Fig. [5b have an enhancement l/m q from the pion propagator at 
zero momentum. The Adler consistency condition guarantees that the 1 / m q enhancement 
in the pole diagrams is canceled, but it does not protect against the finite piece we are 
interested in. 
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Figure 5: Examples of contact a) and pole b) diagrams for the induced 7r~NN vertex. The 
pole diagrams may contain non-singular terms as well and these are included in the contact part 
of the amplitude. The light shaded blob represents the pion amplitudes off the nucleon, the 
solid block shows the insertion of the CP-odd operator 0^~\ The dashed-dotted line is the 7r° 
propagator with an infinitesimal momentum. 



The conventional PCAC vertex derived from the axial charge commutator, the first 
term in Eq. flA.lj) . is just the above contact vertex. To determine the extra finite part we 



take ^=ui^z ) u — di^<~ ) d. Using the operator identity 

1 

ui^u-di^d = — d ll J fl5 

Tflq 
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we have 

(Nir\(m-f 5 u-di^d)(0)\N) = i^(N7r\J* 5 (0)\N) = at g M ^0. (A.8) 

Tflq 

This equation is valid for arbitrary (even large!) nonzero m q and arbitrary ir~ momentum. 

We can examine it to the first two orders in m q . Since there is a pion propagator pole, 
the leading constraint is the vanishing of the l/m q piece. This is the Adler consistency 
condition: the pole residue proportional to the ir° emission amplitude at zero momentum 
vanishes. Vanishing of the 0{m° q ) term means the cancellation of the two contributions, 
the pole and the contact amplitudes; no further terms appear due to the Adler condition 
established at the previous step. The latter comes from various regular terms not con- 
taining poles or kinematic singularities and therefore can be calculated simply at m q = 
and p n - —7-0. The standard PCAC commutator applicable for soft tt~ is just this piece. 
The former contribution is the new chirality- violating terms (xm q which spoil the exact 
Adler cancellation of the amplitude in the chiral limit. Eq. ( 1A.8|) fixes it to be exactly 
minus the contact amplitude: 

(0\(m l5 u-di l5 d)(0)\n°) ■ A( - N7T = -(N7i-\(m l5 u-dt l5 d)(0)\N). (A.9) 

This relation is exact in the limit of p^o — > and is similar in spirit to the Tomozawa- 
Weinberg formula [22], yet is simpler and differs since only one pion is soft. It actually 
applies to any hadron state, not only Nti~ . 

Now we can go back to the case of a general 0^~\ It cannot anymore be represented as 
a total derivative, and the matrix element (Nn~\0^\N) does not vanish. It is still given 
by the sum of the contact and the pole diagrams. The former would again be given, for 
soft 7r~, by the Goldberger-Treiman commutator; it should be taken over the unperturbed 
nucleons, since the quark masses are kept nonzero. The latter, the chirally enhanced pole 
diagrams with the strong vertices corrected at order m q depend only on m q but not on 
0^~\ i.e. they are given by QCD proper. The CP-odd operator O* - -* enters them only 
at the tadpole (0|O^ _ ^(0)|7r°), Figs. [5b. Multiplying the tadpole by the strong amplitude 
A(Nit~ — > Ntt°) over from Eq. (IA.9I) we get, for the CP-odd part, 

(Nn-\O^(0))\N) - /n |,- (0|Q( l (0) '!l| 0N " (N7r-\(ui l5 u- di l5 d)(0)\N). (A.10) 

This representation has an advantage of still being valid at arbitrary 7r~ momentum, 
yet only to the leading order in m q (which here means discarding m 9 ///hadr)- It clearly 
conforms to Eq. flA.7j) . 

The explicit form of the N — > Ntc amplitudes at arbitrary pion momentum is not 
known, therefore to have a concrete expression we finally should assume p w - ^C/ihadr- Using 
the generalized Goldberger-Treiman relation for the conventional unperturbed nucleon 
states (m q is finite now) we expectedly arrive at Eq. (lA.lj) . Thus, we have traced how the 
chiral pole resummation generates the subtraction term exactly in the way anticipated in 
our original simple derivation. 

The additional general observation is useful in view of the left-handed structure of 
the weak currents in the SM. Namely, for any operator not containing right-handed u- 
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or rf-quark fields the induced irNN coupling at zero pion momentum must vanish in the 
chiral limit m u ^— >0. In the cases where the commutator with the axial charge does not 
vanish explicitly this implies the vanishing of the corresponding zero-momentum matrix 
element. This applies to any on-shell amplitude off the hadrons, not only to the nucleon 
vertex, and is a counterpart of the Adler consistency condition. 

The reason is that at m Uj d = the theory is invariant under the isotriplet right-handed 
chiral transformation 

q(x)^e i ^ 5)T3 q(x); 

as long as C w is free from ur and d R this symmetry persists in the full theory includ- 
ing the weak interaction. Since it is spontaneously broken by the conventional strong 
dynamics, there is an exactly massless (at m q = 0) Goldstone boson, ir°, associated with 
the corresponding exactly conserved Noether current. This current is evidently a sum of 
the usually considered axial current and of a flavor- diagonal vector current. Likewise, as 
in the conventional axial-current case, the corresponding Goldstone vertex off the exact 
eigenstates vanishes at zero momentum. 

The formal derivation is straightforward if one considers, instead of the conventional 
axial current, the corresponding left-handed current. Its divergence, the analogue of 
Eq. (1A.4|) vanishes at m q = by virtue of the exact quark field equations of motion, and 



the generalized Goldberger-Treiman relation Eq. (I A. 51) says that c(0) =0, cf. Eq. (IA.6j) . 



The existence itself of the exact Goldstone boson follows from Eq. flA.5j) considered in the 
limit of small nonvanishing m q with P = 0. 

Therefore, any weak pion amplitude vanishes for small pion momentum in the chiral 
limit unless weak Lagrangian contains right-handed u or d fields. More generally, it may 
only remain finite if there is no combination of vector and non-anomalous axial transfor- 
mation that leaves C w invariant. Furthermore, the invariance can be checked at arbitrary 
normalization scale, and usually it is most evident for the bare operators. As an exam- 
ple, the bare C w in the Standard Model contains only left-handed fields, but Penguins 
[I] induce the operators with the right-handed light quarks in the conventionally consid- 
ered effective renormalized Lagrangian. Nevertheless the zero-momentum pion amplitude 
vanishes in the chiral limit in the Standard Model. 

Unlike the pseudoscalar vertex, the induced CP-violating scalar pion-nucleon vertex 
describes the Lorentz structure in the amplitude that does not vanish at zero momentum 
transfer. Therefore, it must vanish in the chiral limit unless the weak interactions include 
right-handed u or d fields. 

Concluding the brief discussion of the application of the current algebra technique we 
note that the similar methods can be applied, for instance to the usual weak decays, e.g. of 
kaons or hyperons, both parity-conserving and parity-violating. In the parity-conserving 
AS = 1 decays we may subtract from the weak Lagrangian the scalar operator sd with 
an arbitrary coefficient A, rewriting it as d^sj^d) / ' (m s —md). This demonstrates that the 
decay amplitude is not changed regardless of A. For parity- violating transitions we can 
subtract sry 5 <i = d^s^^d) / ' {m s + md). This is also useful in establishing the absence of 
the chiral enhancement in the .fT-decay amplitudes mediated by composite quark bilinears 
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like sia^G^d, and to elucidate other similar cancellations. We do not expand on the 
related applications here. 
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